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We present a calculation of generalized baryon form factors in the framework of three-flavor 
covariant baryon chiral perturbation theory at leading one-loop order, necessary for the calculation 
of the first moments of generahzed parton distribution functions. The formulae we derive can be 
used to guide the chiral extrapolation of three-flavor lattice calculations of the corresponding QCD 
matrix elements. 
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I. INTRODUCTION 



tors 



The concept of generalized parton distributions 
(GPDs) connects several different, seemingly uncon- 
nected hadron structure observables such as form fac- 
tors, angular momentum carried by quarks and gluons, 
moments of parton distribution funtions, transverse spa- 
tial structure, etc.. Thus, GPDs and Mellin moments 
thereof are important quantities for studying hadron 
structure [Hll]. 

Lattice QCD has proven to be able to probe hadron struc- 
ture in the non-perturbative regime quite successfully for 
two dynamical quark flavors [3], but today most simula- 
tions are carried out for three or four dynamical flavors. 
The quark masses in these simulations are typically un- 
physically large as it often seems more efficient to in- 
vest computing time in simulating with larger volumes, 
finer lattices or better fermion actions. Covariant baryon 
chiral perturbation theory (BChPT) [H U] is an effec- 
tive field theory of QCD which supplies extrapolation 
formulae for variable quark masses (pseudoscalar meson 
masses) which are of vital importance to thoroughly an- 
alyze lattice data. 

In this article, we generalize the work presented in ^ 
to three quark flavors. In sect. 2 we shortly review the 
connection between GPDs and the so-called generalized 
form factors. In sect. 3 we present the effective chiral 
Lagrangian needed for a leading one-loop calculation of 
the SU{3) version of in the forward limit and in sect. 
4 we show the results for the nucleon sector. We give a 
short conclusion in sect. 5. 



II. GPDS AND GENERALIZED FORM 
FACTORS 



It is well known that the parity-even generalized par- 
ton distribution functions A^) and E'^{x,^, A'^) 
(with Bjorken-x, skewedness ^ and momentum trans- 
fer A) are connected to three generalized form factors 
A| o(A2), B|o(A2) and C|_o(A2) via their first Mellin 
moments [2]. These generalized form factors themselves 
are accessible through calculation of baryon matrix ele- 
ments of totally symmetrized and traceless local opera- 



Here, we have introduced the abbreviations 

A{^B^} = ^ (gat^gfiiy + gp^^Qau - 29afS9tJ.u ] A°' , 



(1) 



(2) 
(3) 



where d represents the space-time dimension. For two 
light quark flavors and assuming isospin symmetry, this 
matrix element can be decomposed into said generalized 
form factors. In the S'C/(3) case, however, one finds five 
different form factors if one imposes no restrictions on 
isospin or baryon content. In this work, we analyze the 
flavor-singlet and the flavor-octet sector, i.e. 



M%,B = {B',s',Y>'\iqt"i{^D^}q\B,s,p 



(4) 



= (S',s',p'|igA^7{M^.}'?|S,s,p), (5) 

where the A' denote the Gell-Mann matrices and B{B') 
labels the incoming (outgoing) baryon from the lowest- 
lying baryon octet. The decomposition of these matrix 
elements yields the following result: 
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Pl^A,.} + £;b',s(A2)7{^A^} 



u{p). 
(6) 



Here, we have introduced another standard momentum 
variable p = {p' + p)/2. Moreover, m = (ms -I- mB')/^, . 
Taking C- and T-symmetry into account, we find that 
A*'", B'^'^ and C'"" are Hermitian 8 x 8- matrices, whereas 
D'^'^ and E'^''" are anti-Hermitian. These matrices are di- 
rectly accessible via SU{3) BChPT. 
As also mentioned in [5], in the forward limit, the gener- 
alized form factor q (0) is linked to the first moment of 
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the parton distribution functions (PDFs) q{x) and q{x) 
via the relation 



{x^ = / dx X [q{x) + q{x)] 



(7) 



and thus, certain A^,^(0) or hnear combinations thereof 
are connected to a hnear combination of these (a;)^, e.g. 
for the isovector moment, we find 

A;'^'{0) ^liix)^ - (xy) ^l{x)u-,. (8) 



III. EFFECTIVE LAGRANGIANS 

Chiral Perturbation Theory provides low-energy ex- 
pansions of QCD Green functions in terms of a small 
parameter p (small compared to a typical hadronic scale 
of ~1 GeV), where p can stand for meson four-momenta, 
baryon three-momenta or meson masses. The well-known 
leading order SU{3) Lagrangian [Tj in the one-baryon- 
sector reads 
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mo{BB) 
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+ -^{BrTo{u^.,B}) + -{Brj,[u,„B]), 
where B denotes the baryon octet 
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mo represents the baryon mass in the chiral limit and 
D and F are the two axial-vector coupling constants. 
Furthermore, we collect the pseudoscalar fields in a 3 x 3 
unitary matrix u which is defined as w = exp{i$/i^o} j 
where 
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contains the pseudoscalar octet of (pseudo-) Goldstone 
bosons associated with the spontaneously broken approx- 
imate chiral symmetry of QCD, and Fq is the perti- 
nent meson decay constant in the chiral limit. Z?^ is 
the appropriate covariant derivative and is defined as 
D^B = d^B + \r^^B] , with the chiral connection (we 
set external vector or axial-vector fields to zero here) 



(12) 



The operator is defined 
again without external vector or axial-vector fields, and 
transforms in the same way under chiral rotations as the 
matter-field B [71[H]. 



As proposed in [B], we now extend this Lagrangian to 
the interaction between the baryon octet and symmet- 
ric, traceless external flavor-singlet and flavor-octet 
tensor fields v)^^^} of definite parity. For the construction 
of the Lagrangian, we utilize the tensor structures 



(13) 
(14) 
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We do not assign a chiral power to these tensor structures 
and thus our leading order Lagrangian is of zeroth order. 
We can now write down all terms that are allowed by 
the symmetry properties of QCD and Lorentz invariance, 
which leads to the following result: 



(0) 



an 



MB,t 



{Bzr{V+,D^B})+h.c. 



ap 



Bir[V+,D-B])+\,.c. 



+ ^{mri,{V-,,D'^B})+\,.c. (15) 
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(B*7^75[K".,i?''S]>+li.c. 



+ ^{BirV°,D"B)+\i.c. 

For the construction of the C'(p^) Lagrangian we need 
additional building blocks, namely the chiral symmetry 
breaking term x which contains the quark mass matrix 
M: 



= diag(m„,md,TOs) , (16) 



X+ = u^xu^ + ux^u, X+=X+- l/3(x+)- (17) 

Using these definitions we can construct the second order 
Lagrangian needed for our leading-loop calculation: 

^MB,t = h{Bir{V+,D^B}){x+)+\^.c. 

+ t2{Bir[V+,D''B]){x+)+\^.c. 

+ u{mr{{v+, X+}, D'-B}) + h.c. 

+ U{Bir[{V+,x+},D^B])+h.c. 

+ h{Bir{K,,X+lD''B})+h.c. 

+ tQ{Bir[[V+,x+lD''B])+h.c. 

+ t,{Bir{V+, {x+, D'^B}}) + h.c. (18) 

+ h{Bir[V+,{x+.D- B}]) + h.c. 

+ h{BirD''B){V+x+)+\^.c. 

+ ho{BirV+){x+D''B)+h.c. 

+ tn{Btj'^{V^,.,D''B}){x+) +h.c. 

+ tu{Btr{{V;^,,X+},D''B}) +h.c. 

+ t,,{B^J^[{V^^,X+},D''B]) +h.c. 

Note that the coupling tio multiplies a non-hermitian 
strucure and can therefore be complex in general. Also, 
this coupling will give a contribution to the form factor 
E defined in eq. 
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IV. RESULTS 



of the form factors, and in particular for the study of fi- 
nite volume effects, the explicit inclusion of the decuplet 
fields will probably be inevitable in a three-flavor calcu- 
lation. These considerations are still work in progress. 
The standard BChPT power counting formula [5, can be 
written down as 

D^2L + l^ J2{n - + - 1)^MB , (19) 



(a) 



(b) 



(c) 



(d) 



FIG. 1; Feynman diagrams contributing to leading-loop or- 
der: the solid line denotes baryons, the dashed line denotes 
pseudoscalar mesons and the dot represents an external ten- 
sor field from the ©(p") La grangian coupling to the baryons 
and mesons. The tree-level diagrams are not shown here. 



As was mentioned in the introduction, our calculation 
is intended as a direct generalization of [6] to the three- 
flavor case. Therefore, the only fields appearing as dy- 
namical degrees of freedom are the mesons and baryons 
from the ground-state octets of flavor-S'[/(3). The effects 
of all other fields are (at least to the order we are work- 
ing here) encoded in the (so far undetermined) LECs ap- 
pearing in the effective Lagrangians constructed above, in 
analogy to the treatment of the baryon masses in [9l [10] . 
Of course, the problem is more urgent in a three-flavor 
calculation: The threshold energy for e.g. a state, 
which is included in our calculation, is at about 1685 
MeV, while that of ttA is at 1370 MeV, the 7r-Roper 
threshold at roughly 1580 MeV, and so on. Let us con- 
centrate on the decuplet, which would certainly be the 
most important resonant state due to the small — A 
mass splitting. It has been incorporated in many stud- 
ies employing three- flavor HBChPT, see e.g. [TTJ W2\ . 
and an extension of Infrared Regularization to the case 
of dynamical decuplet fields has also been given in the 
meantime |13j . For the observables calculated here, how- 
ever, the inclusion of explicit decuplet fields would not 
render the theory more effective, since again a host of 
new undetermined parameters would have to be intro- 
duced (e.g. yAMi?, 1/AA couplings, where V denotes 
the tensor source field) , of which there are already enough 
to absorb the leading effects of resonances in our present 
framework. It might therefore happen that some of our 
LECs will be enhanced (comparing to natural size esti- 
mates) due to important resonance contributions. On the 
other hand, in order to study the momentum dependence 



where D represents the chiral dimension of a Feynman 
diagram, L denotes the number of loops and ^^mmb 
counts the number of vertices that stem from the me- 
son and meson-baryon Lagrangians of power respec- 
tively. In our case, one has to take into account the fact 
that the lowest order meson-baryon Lagrangian including 
symmetric tensor fields starts at order 0{^) and thus, 
all leading one-loop order diagrams are of ©(p^). These 
topologies are depicted in fig. [T] 

From our calculation, we can extract the flavor-singlet 
and flavor-octet generalized form factor A^^l^ to 0{p^) in 
the forward case for different B' and B. In this work, we 
only show the results for the nucleon-channels, because 
all expressions are rather lengthy and the complete list is 
only of interest for practitioners of such calculations. The 
complete set of results is, however, available electroni- 
cally [14,. We note that some of these matrix elements 
have also been calculated in the framework of partially 
quenched ChPT, see refs. [13 [TS]. Our results are ex- 
act to leading one-loop order 0{p'^) . In our formulae, 
however, we display the full loop functions, which also 
contain terms of higher order in the meson mass (quark 
mass) expansion. We have observed that the truncated 
leading-one-loop results approximate the full loop func- 
tions rather badly for meson masses already above ^ 300 
MeV. In particular, the nonanalytic M'^-contributions to 
the loops are far from negligible. As an example, we 
show the baryon one-loop wave function renormalization 
factor at = Mk = = Msymm in fig. [2| and the 
truncated O(p^) (dashed) and 0{p^) (dotted) approxi- 
mations to it. For the other graphs, the situation is very 
similar. Therefore, to take the numerically important 
higher-order parts into account, we insist to use the full 
loop functions everywhere, in accord with the original 
proposal in Tf . As a side remark, we note that the fact 
that Z ^ 1 for meson masses much above 400 MeV may 
cast doubts on the applicability of BChPT in this regime. 
Projecting our effective Lagrangian onto the SU{2) sec- 
tor, and comparing our results for the one-loop ampli- 
tudes with those derived in [B], we also obtain matching 
relations for the SU{2) parameters g and c'gl to the 

SU{3) parameters au, ap, and . 

The results for the generalized form factors Ag"g{0) 
are presented in eqs. (24)-(27), while the matching rela- 
tions can be found in eqs. (20)- (23). We have left out 



(r) 

the index r of the renormalized couplings tl for bet- 
ter legibility. In the matching relations, one could also 
replace the strange quark mass by the Gell-Mann-Oakes- 
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0.2 0.4 0.6 

A4y,„„ [GeV] 

FIG. 2: One-loop wave function renormalization factor for 
Mt, = Mk = Mr) = Msymm'- fuU rcsult (solid line) and trun- 
cated results (dashed and dotted lines). 



Renner-relation 2BQms = 2M^ — to the order we 
are working in. Note that the singlet form factors do not 
get contributions from loop graphs at leading one-loop 
order. Also, note that A^'^ = A^'^^ , A^'^ = zA^jf , A^:^ = 

_ Av,3 Av.,S _ Av,S 
pp 1 nn pp ■ 



Al^{Q) 



aD + aF + 16Bo"^s(^l + ^2 - ^ (2(^3 + ^4) - ^7 - h)) 



M' 



K 



487r2F2 



apiD'^ - 6DF - 3F^) + ao (~7D^ + 2DF - 



3 {ao (1 + 7D^ - 2DF + 3F^) + ap {l - + 6DF + 3F^)) log (^^) + C(^^i), 
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Bom, (Sill + ti2 - 3ti3) + 0{Ml), 



cl{\ = Too) = ml 4(ti + t2) + -(2(i3 + U) 0{MI) 



C9 = 3^^(6*11 - h2 + 3ii3) + 0{Ml). 



y + ^^4(3^11 + ii2 - 3^13) + Ml (^4tii - ^ + 8ti3 
y + ^Af|.(3iii + ii2 - 3^13) + Ml Uhi - ^ + 8tl3 



(20) 

(21) 
(22) 
(23) 
(24) 
(25) 
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ap +aF 4 
2 3 
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Im{Mk) 



3{h + t2){2Ml + Ml) - 2(2*3 + 2t4 - - U){Ml - M^) 
{D + F){2,{aD + aF){D + F) + 8{AaD + AaF))Ml 

- 3(or, + aF)(4 + {D + Ff)ml 

2{-AaD{D - 3F) - dapiD - Ff + ^ao^D - F){D + 3F) + 3{D + F)AaF)M^ 

- 3(3aF - 6aF{D - Ff + aoi^ + 2{D - F){D + 3f )))to^ 



Imb(Mt,) 
Imb{Mk] 



{D + F)Ml [-16(Aaz5 + AaF)(4mg - M^) + 3(ac + aF){D + F){-Srnl + bMl 



(36F2m2) 



Mi 



2{D{AaD - 3Aair) - 3F(AaD + AaF)){^ml - M%) 
+ 3(D - F){{aD - 2,aF)D + 3(ar, + aF)F){-8ml + 5M|.) 



+ 



(48F2m§) 
Jmbb(0,M^) 
{l&F^mD 

Imbb{0,Mk) 
{l2F^ml) 

lMBBiO,M^] 

{48FSml) 

{D + F)M* 
(23047r2i^o2^^) 
Ml, 



[ao + aF){D + F)2M2(8m^ - l2mlMl + 3M^) 
(D - F){{aD - SaF)D + 3{aD + aF)-F)M|.(8m^ - 12m^M| + 3M|.) 



(ac + air)(i? - 3FfM^{8mt^ - 12mgM^ + 3M^) 

9(az, + of)^ + F){2ml - Ml) + 4(Aar, + AoF)(6m^ - M^) 
MaF{D - F)2(2mg - M^) - 18aD{.D^ + 2DF - 2,F^){2ml - M^) 
+ {D{AaD - iAap) - 3F{AaD + AaF))(6mg - M]^) 



+ 



(34567r2F2m^ 

{an + aF){D - 3F)^2mlM^ - M^) 



(7687r2Fo2m^) 
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N 



au — Zap 4 
2^/3 3\/3 

Im{Mk) 



Ziti - it2){2Ml + Ml) + 2(10^3 - 6t4 + iy - Sig + 6f9)(M|. - Ml) 



+ - iap){D + F)'{-ml + Ml) 



(UV^Fgrnl) 



2{5DAaD - AaDD{D - 2,F) ~ SAaoF - SDAap + 9FAaF)M]^ 
+ {-27aF + ac(9 + 8D{D - 3F)))nil 



Im {M„) _ ^^^^^^ _ ^pY^„,2 „ ^^2) 



(24V3F2m2) 



lmb 



Ml 



{IQF^ml 
Imb{Mk) 
ieVSF^ml) 

iASVSF^ml) 
Imbb{0,M^) ,- 
{16F,^m.l) ""'^ 

iMBBiO, Mk 



V3{aD -3aF)iD + F)^Ml{8ml - 5MI) 



(3x/3F2m§) 

/mbb(0,M^) 



{5DAaD - 'SAaoF - SDAap + 9-FAai.)(4mg - M|.) 
+ 2aDD{D - ?,F){-%ml + 5M|.) 

as - 'iaF){D + i^)2Af2(8TO;^ - 12to2m2 ^ 3^4^ 



(flD - 'iaF){D - 3F)2Af2(8m^ - 12mgA/^2 + SM,'^) 



(48v^F2m2) 

x/3(az3 - iaF){D + f )^(2mgAr4 - M^) 
(2567r2F2m4) 

- {bDAao - iAaoF - SDAqf + 9FAaF){6ml - Af|.) 

+ 24:aDD{D - 3F){2ml - M^) 
{ao - 3aF){D - iF)^{2mlM^ - M^) 



{llb2V^ii'^F^ml) 



(768\/37r2F2m4) 



(27) 



V. DISCUSSION AND CONCLUSION 

In this work, we have calculated the flavor-singlet and 
flavor-octet generalized form factor A'^^g in the forward 
limit for three dynamical quark flavors. We have pre- 
sented the results for the nucleon-channels including the 
full loop functions, employing the covariant regulariza- 
tion scheme of ^T7] . We have also calculated the matching 
relations between the S'J7(3) parameters and the SU{2) 
parameters introduced in [6]. Together with the infor- 
mation from two-flavor lattice simulations and the cor- 
responding results of [B], these relations can be used as 
constraints on fitting parameters when chiral extrapo- 
lations of three-flavor lattice calculations of the matrix 
elements in eqs. Q, ^ are performed. However, a word 
of caution is in order here. It is well-known that three- 



flavor baryon ChPT is only poorly converging, mostly 
due to the large kaon and eta masses, see e.g. [21 (TB]. It 
is therefore reasonable to expect that a leading-one-loop 
calculation as presented here might not be sufficient for a 
satisfactory description of lattice data near the physical 
regime. 

Despite this known difficulty, there exist some approaches 
to lattice simulations where the chiral extrapolation for- 
mulae derived here should be useful. The strategy fol- 
lowed in [ini HO] is to start from the S'C/(3)-symmetric 
point where m„ = — ms and to approach the physi- 
cal point keeping the average (singlet) quark mass fixed. 
For realistic values of the quark masses, this corresponds 
to octet meson masses of M^y^^ = i(2Af|- -I- Ml) w 
(411 MeV)^. So, in the vicinity of the symmetric point, 
the kaon and eta masses are considerably smaller than 
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in the real world. As one can easily see from eq. (18), 



only the constants ti, t2, in are relevant at the symmetric 
point, while the other couplings only parametrize devi- 
ations from the SU{3) symmetry limit. In the case of 
baryon masses, it has been argued in [21j that an accu- 
rate extrapolation to the physical point can be obtained 
already employing an expansion linear in the symmetry 
breaking due to the quark masses, which amounts to a 
linear extrapolation in M'^y^^ — . This seems to work 
reasonably well and leads to so-called fan plots for the 
baryon masses, shown in [3T]. In the case of GPDs, how- 
ever, one should expect that it is more important to take 
the proper chirally nonanalytic behaviour for small pion 
masses into account, because already the leading correc- 
tion of 0{p^) to these form factors features chiral loga- 
rithms. This is not the case for baryon masses, where the 
leading quark mass correction is just given by mass inser- 
tions from local counterterms in the chiral Lagrangian. 
In the following, we would like to argue that the men- 
tioned type of lattice simulations is well suited for an 
application of our extrapolation formulae, presented in 
the last section. To get a rough estimate of the size of 
possible higher-order corrections to these results for this 
case, we consider the contributions from the graphs in 
figs. 1(a) and l(b)[ which turn out to start only at next- 
to- leading loop order 0{p^). They are proportional to 
the couplings AaD,F, of which only the order of magni- 
tude is (roughly) known, see table 2 of !6j for the value 
of the corresponding SU{2) parameter Aajg. We vary 
these couplings in the range —0.3 < l^ajj,/S.aF < 0.5, 
with all other couplings fixed at typical values consistent 
with the analysis in the SU{2) sector and earlier work 
on BChPT {ao ^ ap ^ 0.1, mo = 1 GeV, D = 0.8, 
F = 0.5 and Fq = 0.09 GeV) and aU ti set to zero. The 
generated band shown in fig. [3] should give a good im- 
pression of the expected size of higher-order contribu- 
tions. In the lower panel of fig. [3) the singlet quark mass 
is kept fixed, while the pion mass is varied, in analogy 
to the "fan plots" shown in [3T]. Regarding the rather 
generous range for the variation of the LECs Aa/j, Aap, 
and earlier experiences with the convergence properties 
of three-flavor chiral expansions, the band in the lower 
panel indicates an improved stability to higher order cor- 
rections, which leaves us optimistic that the presently 
derived chiral extrapolation formulae provide a realistic 
description for the special type of lattice simulations de- 
scribed above. In principle it is straightforward to extend 
the analysis to a next-to-leading one-loop calculation. On 
the other hand, some additional coupling constants will 
appear if one pushes the calculation to higher orders, and 
one would have to specify the individual baryon masses 
in the loop functions, which complicates the calculations 



a lot. Even at the order we are presently working, it will 
require an enormous amount of three-flavor lattice data 

(r) 

to be able to pin down all of our LECst]''. Therefore, we 
defer the extension of the present results to a later com- 
munication. In our opinion, the most promising strategy 
for the moment is to stay with the results as derived here. 
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and to first concentrate on a better understanding of the 
two-flavor sector, e.g. by analyzing higher-order correc- 
tions to the existing calculations '6', '22', ^23| . There, the 
number of new undetermined parameters is much more 
limited, and convergence is a less severe issue than in 
the three-flavor framework. Work along this direction is 
already in progress ^24^ i25j . 
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Appendix A: Loop functions 



In sec. IV we have used the following abbreviations for the loop functions: 

Im{M) = — Ina, 



(Al) 



where a = M/nio ■ Note that we use fi = ttiq for the renormalization scale everywhere, following the original proposal 
of Infrared Regularization by Becher and Leutwyler. Furthermore, 



(M) 



167r2 



(21na- 1) 



arccos — — 



(-1) 



Moreover, the renormalized three-point function in Infrared Regularization, taken at — 0, is given by 



MBB 



1 



327r2TOQ 



2 In a + 1 - 



2a 



: arccos — — 



\/4 - a 

The nucleon wave function renormalization factor, at the one-loop level, is given by 

5D^ - 6DF + 9F^ 



(-1) 



Z 



N 



in^F§) 



" (967r2F2) 



3{D + F)^M^{-3ml + M^) 



(167r2F>3./4-^) 



arccos 



AU\ _ {D~3FfM^{-3ml + M ^) 



{48TT^F§mlJA 



arccos 



Mr, 

2 mo 



{5D'^ - 6DF + 9F2)M|-(-3m2 + M^) 



{2^TT^F^ml.U 



Ml 



arccos 



) 



Mk\ 2,{D + FfMl{-Zml + 2Ml) 



2mo 



m^^F^ml) 



log 



mo 



+ 



{hD^ -QDF + <^F^)Ml{-2,ml + 2Ml)^^^( Mk\ , {D - 'iFf M^(~?>ml + 2M^) 



in^F^ml) 



mo 



{9671^ Fgml) 



log 



mo 



(A2) 



(A3) 
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